7С 


УДК 621.039.55 


А.В. Кочергин, С.С. Пивоварцев 

Научно-исследовательский и проектно-конструкторский институт «Искра», 
г. Луганск, Украина 

Кау@пе*.1е ма 


Нейронная сеть для идентификации нуклидов 
по гамма-спектру 


В статье представлен метод идентификации радионуклидов по реконструированному физическому 
спектру гамма-излучения с использованием трехслойной нейронной сети прямого распространения. 
Описана структура сети, способы формирования обучающей выборки и результаты экспериментальной 
проверки метода. 


При оперативном контроле и анализе радиационной обстановки непосредствен- 
но на местности в режиме реального времени не достаточно измерить дозовые 
характеристики поля. Важно также определить потоковые характеристики поля, 
спектральный состав излучения и радионуклидный состав загрязнения с оценкой 
активности идентифицированных изотопов. 

Радиометрическая обработка полученных при экспериментах спектров гамма- 
излучения, как правило, заключается в извлечении информации о положении и интен- 
сивности спектральных линий и их идентификации по энергиям с радионуклидами 
каталога. Эта информация чаще всего может быть получена при изучении наиболее 
информативной части спектров — пиков полного поглощения (ШП). Положение ППП в 
спектре дает информацию об энергии данной линии, а его площадь — информацию об 
интенсивности. В последнее время наблюдается тенденция к применению для решения 
подобных задач нейроподобных структур [1]. В рамках данной работы представлена 
нейронная сеть для идентификации нуклидов по гамма-спектру, реализованная в 
портативном полевом приборе для обнаружения незаконного или непреднамеренного 
перемещения ядерных или радиоактивных материалов. Базовые требования к прибо- 
рам этого класса определены международным стандартом ТЕС 62327 [2]. 


Постановка задачи 


В общем виде задача выглядит следующим образом. 

Основываясь на данных аппаратурного спектра, необходимо определить качест- 
венный радионуклидный состав исследуемого радиоактивного вещества и приблизи- 
тельные массовые доли входящих в вещество нуклидов. 

Данная задача в терминологии нейронных сетей относится к классу классифика- 
ционных (райеги с1азсаноп). Для решения этого класса задач с успехом применяются 
многослойные полносвязные сети прямого распространения. 

Стандартом определен следующий список основных радионуклидов, которые 
подлежат идентификации. Ядерные материалы: ?°°0, 250, ??°Ри, "Ри. 

Медицинские радионуклиды: 13Е, °'ба, Те, ИЧа, 1, Тодте ( о" у 
201]. 133Хе, природные радионуклиды (МОЕКМ): К, 2*Ва, ТВ, 330; технические 
радионуклиды: 16, 5: 33Ва, Е: 92, Ва, 21 Ат, О РЫ: 
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Лишь незначительная часть радионуклидов из указанного списка имеют относи- 
тельно простой спектр гамма-линий. Условия и временные ограничения процедуры 
измерения и идентификации приводят к тому, что в большинстве случаев резуль- 
тирующий гамма-спектр имеет сложный характер со значительным количеством 
мультиплетов, тяжелых для идентификации при реальных уровнях скорости счета, 
не намного превышающих фоновую загрузку детектора. На рис. 1 представлен 
суммарный линейчатый спектр энергий гамма-линий с учетом вероятности выхода 
для всех изотопов предложенной библиотеки радионуклидов. 


о 0.5 у] = 7. 25 з 
Е Ме\ 


Рисунок 1 — Линейчатый гамма-спектр изотопов библиотеки нуклидов 


Учитывая этот факт и относительное энергетическое разрешение возможных 
по условиям эксплуатации детекторов, очевидно, что число одновременно иденти- 
фицированных изотопов не может быть велико. Основываясь на требованиях 
стандарта, а также учитывая тот факт, что в фоновом спектре всегда присутствует 
определенный набор природных радионуклидов из списка МОКМ, было принято 
решение ограничить число одновременно идентифицируемых нуклидов пятью. 

Для исследуемых образцов, которые представляют собой смесь элементарных 
источников, получаемый спектр 5 может быть аппроксимирован линейной суперпо- 
зицией его составляющих 5;: 5 = а .5,. Отсюда очевидно, что система идентифика- 


1 


ции должна иметь линейную зависимость выхода от входа. И хотя обычно используют 
свойство нелинейности ИНС, преимущества такого подхода неоспоримы ввиду 
простоты, скорости и автоматизации принятия решения. 

Попытка построить и обучить одну общую нейронную сеть для одновременной 
однопроходной идентификации любого из нуклидов из определенной выше библио- 
теки столкнулась со значительными трудностями и не привела к удовлетворяющему 
нас результату. Путем последовательных приближений было принято решение перейти к 
многопроходной процедуре идентификации, когда для поиска каждого нуклида из 
библиотеки строится и тренируется отдельная нейронная сеть. 


Методика обучения сети 


Входными данными для нейронной сети являлся реконструированный физический 
256 канальный спектр гамма-излучения в точке измерения. Под восстановленным 
физическим спектром понимается энергетическое распределение плотности потока 
гамма-квантов в данной точке, оцененное по аппаратному спектру. 
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Задача восстановления физического спектра в точке измерения по измеренному 
аппаратному спектру может быть сведена к решению матричного уравнения (1) 


А:и=Х, (1) 


где А — представляет собой матрицу аппаратной функции измерительного тракта, /— 
измеренный аппаратный спектр, и — физический гамма-спектр в точке измерения. 
Измеренный спектр /, как и аппаратная функция А, известны нам с некоторой 
погрешностью. При этом возникают трудности, связанные с некорректностью задачи 
по Адамару. Для решения уравнения (1) был применен итерационный метод Скоффилда 
Голда [3], суть которого заключается в следующем. 

Уравнение (1) преобразуется к виду 


В-фЕЙ, (2) 
где В=А'.А; Т=А'.Х. 


На каждом шаге п итерационной процедуры вычисляется вектор Г” =В-и" и 


п 


И. 
строится диагональная матрица с элементами 4, =—. 
› у" 
1 
Следующее приближение определяется по формуле 
и" =Р".И’. (3) 


Данная процедура позволяет существенно повысить разрешающую способ- 
ность и точность предлагаемого метода. 

Процедура реконструкции проводилась по модельному аппаратурному спектру. 
Применение модельных, а не реальных аппаратурных спектров определялось прак- 
тической невозможностью экспериментального получения спектров всех возможных 
комбинаций нуклидов и их относительных активностей. В связи с этим в основу 
получения входных данных были положены две методики. 

1. Реализованная на базе программного комплекса СЕАМТ 4, разработанного в 
Европейском центре ядерных исследований (СЕКМ, Женева), программная модель 
для имитации взаимодействия частиц и излучения с веществом. Разработанная нами 
на базе СЕАМТ программная модель взаимодействия поля излучения с произволь- 
ным детектором, позволяет оперативно изменять как количество и форму источни- 
ков излучения и их радионуклидный состав, так и параметры детекторов и среды 
измерения. 

2. Спектры, полученные умножением вектора ожидаемого линейчатого спектра 
гамма-излучения в точке измерения с матрицей аппаратной функции спектрометри- 
ческого тракта размером 1011х1011 элементов сгенерированной методами Монте- 
Карло по схеме индивидуальных соударений программой «ТгапзЁег» разработки 
НИПКИ «Искра». 

Кроме этого в модельный спектр добавлялись фон, статистический шум и слу- 
чайный дрейф по шкале канал — энергия в диапазоне + 1,5 %. 

Предложенная методика позволяет создать базу данных спектров гамма-излучения, 
охватывающую все типовые ситуации обнаружения радионуклидов, изложенные в 
требованиях действующих стандартов и дополняющую имеющиеся эксперименталь- 
ные данные. 

Проверка предложенных алгоритмов проводилась как на модельных спектрах, 
не входящих в обучающую выборку, так и на реальных спектрах, полученных на 
спектрометрической аппаратуре с детектором на базе сцинтилляционного кристалла 
МакТу ©25х30 мм + ФЭУ и С$(Т] ©40 х 40 мм + ФЭУ. 
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ТК 


Результаты работы сети 


Построение и тренировка сети осуществлялась с помощью средств Меига] МебмогК$ 
Тоофох МАТГАВ 7.6. 
Структура сети представлена на рис. 2. 


256 256 10 1 


Рисунок 2 — Структура нейронной сети 


Для выбора числа нейронов в скрытом слое использовалась формула для оценки 
необходимого числа синаптических весов №, в многослойной сети с сигмоидальными 
передаточными функциями 


г <М№ а :(м +М№ +1)+№ 
1+8, (№, “ом, о й 


где №, — размерность выходного сигнала, №, — число элементов обучающей выборки, 
№. — размерность входного сигнала. 

Оценив необходимое число весов, можно рассчитать число (№ нейронов в 
скрытых слоях 


№ 


ГМ, М, 


В любом случае расчетная величина М является оценочным параметром, имею- 
щим существенный диапазон варьирования. 
Для обучения сети был применен метод сопряженных градиентов с адаптивной 
корректировкой длины шага (функция га!1$се). 
На рис. 3 и в табл. 1 представлены результаты работы сети для тестовой комби- 
нации нуклидов и их относительных активностей в спектре. 
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Исходный и фоновый спектры Восстановленный спектр 


Рисунок 3 — Спектры гамма-излучения 
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Таблица 1 
Изотоп Доля в Превышение над Результат работы 

спектре фоном, % сети 

Ва-133 0 0,0066 

Со-60 0,15 0,9813 

С5$-137 0,65 75 0,9760 

ТЬ-228 0,2 0,9983 

0-238 0 0,0197 


Оценка результатов 


Резюмируя проведенные исследования, можно сделать следующие выводы. 

Метод показал свою безусловную перспективность, в особенности в применении 
к сцинтилляционной гамма-спектрометрии, для которой характерно существенное 
количество мультиплетов, затрудненных к локализации и идентификации из-за естест- 
венных ограничений сцинтилляционных детекторов, связанных с относительно плохим 
энергетическим разрешением. Несомненным преимуществом метода является тот факт, 
что процедура идентификации выполняется за один проход, что безусловно важно, 
принимая во внимание ограниченные возможности вычислительного ядра полевого 
портативного прибора. 

Безусловно, метод достаточно трудоемок на предварительном этапе и весьма 
чувствителен к качеству и количеству информации в обучающей выборке. Подго- 
товка данных требует точного знания технических характеристик и особенностей 
спектрометрического тракта и априорной информации о предполагаемых условиях 
применения разрабатываемого прибора. 

Обученная нами нейронная сеть показала отличные (близкие к 100 %) резуль- 
таты для идентификации единичных нуклидов с простым линейчатым спектром. В то 
же время в случаях сложных спектров либо при идентификации нуклидов, обладающих 
сложным линейчатым спектром (таких как ТВ-232 либо 0-238), результаты неодно- 
значны. Положительные результаты получены для не более чем трех нуклидов, 
одновременно присутствующих в спектре и не во всем диапазоне исследуемых 
относительных активностей. 

Безусловно, метод нуждается в дальнейшем развитии. Перспективным было бы 
ограничение входных данных обучающей выборки теми каналами, в которых гипо- 
тетически могут находиться пики полного поглощения искомого нуклида. 
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